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1 Introduction 



Since the conception by Einstein of General Relativity Theory, a great many efforts have 
been devoted by many scientists to the construction of a consistent quantum theory of 
gravity. These efforts can de divided into two main approaches: 

(a) perturbative, in which one splits the metric into a background (kinematical) part 
and a dynamical one: g^^ = rj^y + h^y and tries to quantize h^y. The only conclusive 
results existing, are that the theory thus obtained is highly nonrenormalizable 0]. 

(b) non perturbative, in which one tries to keep the twofold role of the metric (kinemat- 
ical and dynamical) intact. A hallmark in this direction is canonical quantization. 

In trying to implement this scheme for gravity, one faces the problem of quantizing a 
constrained system. The main steps one has to follow are: 

(i) define the basic operators 'g^^p and t^^'^ and the canonical commutation relation 
they satisfy. 

(ii) define quantum operators whose classical counterparts are the constraint func- 
tions H^. 

(iii) define the quantum states ^[g] as the common null eigenvector of H^, i.e. these 
satisfying H^^[g\ = 0. (As a consequence, one has to check that form a closed 
algebra under the basic CCR.) 

(iv) find the states and define the inner product in the space of these states. 

It is fair to say that the full program has not yet been carried out, although partial steps 
have been made 0. 

In the absence of a full solution to the problem, people have turned to what is 
generally known as quantum cosmology. This is an approximation to quantum gravity 
in which one freezes out all but a finite number of degrees of freedom, and quantizes the 
rest. In this way one is left with a much more manageable problem that is essentially 
quantum mechanics with constraints. Over the years, many models have appeared in the 
literature 0. In most of them, the minisuperspace is flat and the gravitational field is 
represented by no more than three degrees of freedom (generically the three scale factors 
of some anisotropic Bianchi Type model |^). 

In order for the article to be as self consistent as possible, we include in section 2, a 
short introduction to the theory of constrained systems and in section 3, the Kantowski- 
Sachs model as an interdisciplinary example. In section 4, we present the quantization 
of the most general Type II, VI(Class A) & VII(Class A) Vacuum Bianchi Cosmologies. 
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2 Elements of Constrained Dynamics 



2.1 Introduction 

In these short notes, we present the elements of the general methods and some techniques 
of the Constrained Dynamics. It is about a powerful mathematical theory (a method, 
more or less) -primarily developed by P. A. M. Dirac. The scope of it, is to describe 
singular (the definition is to be presented at the next section) physical systems, using a 
generalization of the Hamiltonian or the Lagrangian formalism. This theory, is applicable 
both for discrete (i.e. finite degrees of freedom) and continua (i.e. infinite degrees of 
freedom) systems. 

For the sake of simplicity, the Hamiltonian point view of a physical system is adopted, 
and the discussion will be restricted on discrete systems. A basic bibliography, at which 
the interested reader is strongly suggested to consult, is quoted at the end of these notes. 
Also, the treatment follows reference 0. 



2.2 The Hamiltonian Approach 

Suppose a discrete physical system, whose action integral is: 

A = j Ldt (2.2.2.1) 

The dynamical coordinates, are denoted by g*, with i G [1, . . . , A^] The Lagrangian is a 
function of the coordinates and the velocities, i.e. L = L{q\q^). 

A note is pertinent at this point. If one demands the action integral ( 2.2.2.1| ), to be 
scalar under general coordinate transformations (G.C.T.), then he can be sure that 
the content of the theory to be deduced, will be relativistically covariant even though 
the form of the deduced equations will not be manifestly covariant, on account of the 
appearance of one particular time in a dominant place in the theory (i.e. the time 
variable t occurring already, as soon as one introduces the generalized velocities, in 
consequently the Lagrangian, and finally the Lagrange transformation, in order to pass 
from the Lagrangian, to the Hamiltonian). 

Variation of the action integral, gives the Euler-Lagrange equations of motion: 

d fdL\ dL 



In order to go over to the Hamiltonian formalism, the momentum variables Pi, are 
introduced through: 

dL 

P^=gTr, V i (2.2.2.3) 

In the usual dynamical theories, a very restricting assumption is made; that all momenta 
are independent functions of the velocities, or -in view of the inverse map theorem for a 



3 



function of many variables- that the following (Hessian) determinant: 

\HiA = \—^\ (2.2.2.4) 

is not zero in the whole domain of its definition. If this is the case, then the theorem 
guarantees the validity of the assumption, permits to use the Legendre transformation, 
and the corresponding physical system is called Regular. If this is not the case (i.e. some 
momenta, are not independent functions od the velocities), then there must exist some 
(say M) independent relations of the type: 

0„(g,p) = O, me[l,...,M] (2.2.2.5) 

which are called Primary Constraints. The corresponding physical systems, are charac- 
terized as Singular. 

Variation of the quantity p^g* — L (the Einstein's summation convention is in use), 
results in: 

5 {p,ct -L)=... = {6p,) - (^^^ 6q' (2.2.2.6) 
by virtue of (|2.2.2.2|) . One can see that this variation, involves variations of the g's and 



the p's. So, the quantity under discussion does not involve variation of the velocities and 
thus can be expressed in terms of the g's and the p's, only. This is the Hamiltonian. It 
must be laid stress on the fact that the variations, must respect the restrictions ( p.2.2.5| ), 
i.e. to preserve them -if they are considered as conditions (see, e.g. C. Caratheodory, 
"Calculus of Variations and Partial Differential Equations of the First Order", AMS 
Chelsea (1989)). 

Obviously, the Hamiltonian is not uniquely determined for, zero quantities can be 
added to it. This means that the following: 

Ht = H + u^(t)m (2.2.2.7) 

where m^'s are arbitrary coefficients in the phase space (including the time variable), is 
a valid Hamiltonian too. Variation of ( |2.2.2.7| ) results in: 

= M + u'^p- + term that vanishes as (2.2.2.5) 
p. = -|| - m'"!^ - term that vanishes as (|2.2.2.5|) v • • • ; 

These are the Hamiltonian equations of motion for the system under consideration. This 
scheme, reflects the previous observation about variations under which, conditions must 
be preserved. 

In order to proceed, a generalization of the Poisson Brackets is needed to be intro- 
duced. This is done as follows: 

If /, h are quantities on a space, endowed with a linear map { , } such that: 

{/, g} + {g^ /} = Antisymmetry 

{/ + g:h} = {/, h] + {g, h} Linearity (o o o q-\ 

{fg,h} = f{g,h} + {f,h}g Product Law y^.^.^.^) 

{/, {g, h}} + {g, {h, /}} + {h, {/, g}} = Jacobi Identity 
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If the space is the phase space, then these Generahzed Poisson Brackets, are reduced to 
the usual ones: 

{f^d} dp^ dq^ dpi ^ 

otherwise are subject to the previous algebra -only. 

For a dynamical variable -say g, one can find -with the usage of: 

and of ( ^.2.2.8| ), as well as the generalized Poisson Bracket Algebra ( |2.2.2.9D : 

9-{9,Ht} (2.2.2.12) 

The symbol ~ is the Weak Equality symbol and stands for the following rule (deduced 
from the thorough analysis of the previous procedure): 

A constraint, must not be used before all the Generalized Poisson Brackets, are calcu- 
lated formally (i.e. only with the usage of the algebra (^.2.2.f^ and the usual definition 
( \2.2.2.1C\ ) -when the last is applicable). This rule, is encoded as: 

0^(g,p)^O, mG[l,...,M] (2.2.2.13) 

In the previous procedure, the position of that rule, refiects the need to manipulate the 
■u™'s, which may depend on t only -since they are unknown coefficients, the definition 
( p. 2. 2. 10]) can not be used. 

If the dynamical variable g is any one of the constraints, then (|2.2.2.5| ) declare the 
preservation of zero. Thus, consistency conditions, are deduced: 

{(l)m',H} + u"'{(l)m',<Pm} = (2.2.2.14) 

There are three possibilities: 

CCi Relations ( |2.2.2.14|) lead to identities -maybe, with the help of ( p.2.2.5| ). 

CC2 Relations ( p.2.2.14[ ) lead to equations independent of the u's. These turn to be 
constraints also. They called Secondary, but must be treated on the same footing 
as the primary ones. 

CC3 Relations ( |2.2.2.14[ ) impose conditions on the u's. 

The above procedure must be applied for all the secondary constraints. Again, the 
possible cases will be the previous three. The new constraints which may turn up are 
called secondary too. The procedure is applied for once more and so on. At the end, 
one will have a number of constraints (primary plus secondary) -say J'- and a number 
of conditions on the m's. A detailed analysis of the set of these conditions, shows that: 

= U^{q,p)+V''{t)Vr{q,p) (2.2.2.15) 



5 



where V^{q,p) are the a (in number) independent solutions of the homogeneous systems: 

vr(g,p){0m",0™} = o 

The functions V"(t) are related to the gauge freedom of the physical system. 

Some terminology is needed at this point. 
A dynamical variable R, is said to be First Class, if it has zero Poisson Bracket, with all 
the constraints: 

{R,^n} = 0, ne[l,...,J] (2.2.2.16) 

where is the total number of constraints -i.e. primary plus all the secondary ones. 
It is sufficient for these conditions, to hold weakly -since, by definition, the 0's are the 
only independent quantities that vanish weakly. Otherwise, the variable R, is said to be 
Second Class. If R, is First Class, then the quantity {R, 0^} is strongly equal to some 
linear combination of the (p^s. The following relative theorem (with a trivial proof) holds: 
"The Poisson Bracket of two First Class quantities, is also First Class". 

Using the result ( |2.2.2.15| ) the Hamiltonian ( |2.2.2.7[ ), which is called Total Hamilto- 



nian, is written: 

Ht = H + U^<Pm + V^V^'K = H' + V''<Pa (2.2.2.17) 

with obvious associations. It can be proved that H' and 0a, are first class quantities. 
With this splitting and the relation (|2.2.2.12|) for a dynamical variable g, it can be 
deduced that: 

The First Class Primary Constraints 4>a, are the generating functions (i.e. the quantities 
{g^(pa}) of infinitesimal Contact Transformations; i.e. of transformations which lead to 
changes in the q 's and the p 's, but they do not affect the physical state of the system. 

Successive application of two contact transformations generated by two given First 
Class Primary Constraints and taking into account the order, leads -for the sake of 
consistency- to a new generating function: {g, {(pa, 4'a'}} ■ Thus one can see that First 
Class Secondary Constraints, which may turn up from {(j)a,4>a'}, can also serve as gen- 
erating functions of infinitesimal Contact Transformations. Possibly, another way to 
produce First Class Secondary Constraints, is the First Class quantity {H',(j)a}- Since 
no one has found an example of a First Class Secondary Constraint, which affects the 
physical state when used as generating function, the conclusion is that all First Class 
quantities, are generating functions of infinitesimal Contact Transformations. Thus, the 
total Hamiltonian should be replaced by the Extended Hamiltonian He, defined as: 

He = Hr+Wcpa" (2.2.2.18) 

where the 0a" 's are those First Class Secondary Constraints, which are not already 
included in Ht- Finally, the equation of motion for a dynamical variable g ( |2.2.2.12D is 
altered: 

g-{g,HE] (2.2.2.19) 
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2.3 Quantization of Constrained Systems 
2.3.1 No Second Class Constraints are Present 

The quantization of a classical physical system, whose Lagrangian, gives first class con- 
straints only, is made in three steps: 

5*1 The dynamical coordinates g's and momenta p's, are turned into Hermitian Oper- 
ators q 's and p 's, satisfying the basic commutative algebra: [q^,Pj] = iSj. 

S2 A kind of a Schrodinger equation, is set up. 

5*3 Any dynamical function, become Hermitian Operator -provided that the ordering 
problem is somehow solved. 

Obviously, the constraints -being functions on the phase space- are subject to the 5*3 
rule. Dirac, proposed that when the constraints are turned into operators, they must 
annihilate the wave function \1/: 

0i^ = O, V i (2.2.3.1) 

Successive application of two such given conditions and taking into account the order, 
for sake of consistency, results in: 

[0i,0j]^ = O (2.2.3.2) 

In order for operational conditions (|2.2.3.2| ) not to give new ones on one demands: 

=Cj,0fc (2.2.3.3) 

If it is possible for such an algebra to be deduced, then no new operational conditions on 
\l/ are found and the system is consistent. If this is not the case, the new conditions must 
be taken into account and along with the initial ones, must give closed algebra, otherwise 
the procedure must be continued until a closed algebra is found. The discussion does not 
end here. Consistency between the operational conditions ( |2.2.3.1| ) and the Schrodinger 
equation, is pertinent as well. This lead to: 

[0„i/]^ = O (2.2.3.4) 

and consistency know, reads: 

[0,,^] = (2.2.3.5) 
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2.3.2 Second Class Constraints are Present 

Suppose we have a classical physical system, whose Lagrangian, gives second class con- 
straints. Any set of constraints, can be replaced by a corresponding set of independent 
linear combinations of them. It is thus, in principle, possible to make arrangement such 
that the final set of constraints, contains as much first class constraints as possible. Us- 
ing the remaining -say S in number- second class constraints, the following matrix is 
defined: 

^ij^[Xi,Xj], ii,j)e[l,... ,S] (2.2.3.6) 

where , x's are the remaining (in classical form) second class constraints. A theorem can 
be proved : 

"The determinant of this matrix does not vanish, not even weakly". 

Since the determinant of A is non zero, there is the inverse of this matrix; say A~^. 
Dirac, proposed a new kind of Poisson Bracket, the { , }d' 

s s 

{ , }d = { , } - E E{ ' Xi}A^{x,, } (2.2.3.7) 
i=i j=i 

These Brackets, are antisymmetric, linear in their arguments, obey the product law and 
the Jacobi identity. It holds that: 

{g,HE}D^{9,HE} (2.2.3.8) 

because terms like {xzjHe}, with He being first class, vanish weakly. Thus: 

g!^{g,HE}D (2.2.3.9) 

But: 

{e,xjD = ... = (2.2.3.10) 

if ^ is any of the g's or the p's. Thus, at the classical level, one may put the second class 
constraints equal to zero, before calculating the new Poisson Brackets. That means that: 

Ml The equations x = may be considered as strong equations. 

M2 One, must ignore the corresponding degrees of freedom and 

M3 quantize the rest, according to the general rules, given in the previous section. 
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3 An Interdisciplinary Example: 
The Kantowski-Sachs Model 



The purpose of the present section is to illustrate an application of the Dirac's method 
for constrained systems. The example chosen, is that of Kantowski-Sachs reduced La- 
grangian -i.e. of a vacuum cosmological model; thus the interdisciplinary character of 
the section, emerges. 

Consider, the following Lagrangian: 

^ a{t)h^{t) + 2h{t)a{t)h{t) ^ N{t)a{t) ^^^ ^ ^^ 



2N{t) 2 

where a(t), h{t) and N{t), are the three degrees of freedom (the g's). It must be brought to 
the reader's notice that the Euler-Lagrange equations corresponding to this Lagrangian, 
are tantamount to the Einstein's Field Equations for the Kantowski-Sachs model (de- 
scribed in P]), characterized by the line element: 

ds^ = -N^{t)dt^ + a\t)dr^ + b\t)de^ + b'^{t)sin'^{e)d(f)'^ (3.3.0.2) 

The momenta are: 

dL 

Pn = —■ = 3.3.0.5 

lYiom the third of ( p.3.0.3|) , one can see that there is one primary constraint: 

pjv ~ (3.3.0.6) 

The total Hamiltonian is: 

Ht = H + u{a{t), b{t), N{t),t)pj^ (3.3.0.7) 

where: 

H = pad{t) + pbh{t) -L = N{t)n{t) (3.3.0.8) 

with: 

^^^^ = " 6(ty + W{t) 

The consistency condition (|2.2.2.14| ) applied to: 
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Ai the constraint Pat ~ 0, gives one secondary constraint: 

X = {PN, H} = {pN, N{t)n{t)} = -n{t) ^ (3.3.0.10) 
A straightforward calculation, results in: 

{x,Pn} = (3.3.0.11) 

A2 the previously deduced secondary constraint x ~ 0, gives -by virtue of ( |3.3.0.9| ), 
( p. 3.0.101 ) and ( p.3.0.11| )- no further constraints, since it is identically satisfied 
{CCi case): 

{X, H] + u{a{t), b{t), N{t),t){x,PN} = (3.3.0.12) 



The Poisson Bracket ( |3. 3.0.111) also declares that both p^ and X) fii'st class 
quantities. 

Finally, the equations of motion are: 

d{t) ^ {a{t),HT} (3.3.0.13) 

Pa^{Pa,HT} (3.3.0.14) 

b{t) ^ {b{t),HT} (3.3.0.15) 
Ph^{ph,HT} (3.3.0.16) 
N{t) ^ {N{t),HT} (3.3.0.17) 

Pn^{Pn,Ht} (3.3.0.18) 

The first four equations constitute the usual set of the Euler-Lagrange equations for 
the a{t) and b(t), degrees of freedom. Equation ( p.3.0.17| ), results in the gauge freedom 
related to N(t) since -according to this equation- N(t) = u{a(t),b(t), N(t),t), i.e. an 
arbitrary function of time, while equation ( p. 3.0.181) is trivially satisfied, in view of 
( P.3.0.101) . 
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4 Quantization of the most general Bianchi Type II, 
VI( Class A) & VII( Class A ) Vacuum Cosmologies 



We firstly treat of (see T. Cliristodoulakis, G. O. Papadopoulos, Pliys. Lett. B 501 
(2001) 264-8): 



4.1 The Bianchi Type II Case 

In , we liad considered tlie quantization of an action corresponding to the most general 
Bianchi Type II cosmology, i.e. an action giving Einstein's Field Equations, derived from 
the line element: 

ds'^ = (N^t) - Na{t)N''{t))dt'^ + 2Na{t)af{x)dx'dt + -faf^{t)a^{x)cT^{x)dx'dx(AA.l.l) 
with: 



a"-{x) = a^{x)dx'^ 

dx"^ — x^dx^ 
dx^ 



a^{x 
a'^ix 



a 



\x) = dx^ 



(4.4.1.2) 



da^(x) = IC^y A (jT 



•-"23 



•-^32 



see 
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As is well known [0], the Hamiltonian is H = N{t)HQ + N°'(t)Ha where: 



1 

2' 



(4.4.1.3) 



is the quadratic constraint with: 



7 being the determinant of •jais (the last equality holding only for the Type II case), and: 

Ha = C^aplP,^^' (4.4.1.5) 

are the linear constraints. Note that appearing in the Hamiltonian, is to be identified 
with Nj^. 

The quantities i^o; -f^a, are weakly vanishing i.e. ifo ~ 0, if^ ~ 0. For all class A 
Bianchi Types [C^fi = 0); it can be seen to obey the following first-class algebra: 



{i/o,i^o} = 



(4.4.1.6) 
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which ensures their preservation in time i.e. Hq ^ 0, Ha ~ and estabhshes the 
consistency of the action. 

If we follow Dirac's general proposal [H for quantizing this action, we have to turn 
Hq, Ha, into operators annihilating the wave function 

In the Schrodinger representation: 

^^P^^<^P = _,d (4-4.1.7) 

satisfying the basic Canonical Commutation Relation (CCR) -corresponding to the clas- 
sical ones: 

[%fs,nn = = ^(5^5;; + 6^,6:) (4.4.1.8) 

The quantum version of the 2 independent linear constraints has been used to reduce, 
via the method of characteristics |]12|, the dimension of the initial configuration space 
from 6 (7„/3) to 4 (combinations of 7^/3), i.e. * = ^(^,7,7^2 - 7ii722, 7i27i3 - 711723), 
where q = C^^C^^^la^r^r''- 

According to Kuhaf's and Hajicek's |13| prescription, the "kinetic" part of Hq is to 



be realized as the conformal Laplacian, corresponding to the reduced metric: 

Lap,.^^ = g'' (4.4.1.9) 

where x*, i = 1,2, 3, 4, are the arguments of \E'. The solutions had been presented in 
Note that the first-class algebra satisfied by Hq, Ha, ensures that indeed, all components 
of are functions of the x*'s. The signature of the is (+,+,—,—) signaling the 
existence of gauge degrees of freedom among the x"s. 

Indeed, one can prove that the only gauge invariant quantity which, uniquely 
and irreducibly, characterizes a 3-dimensional geometry admitting the Type II symmetry 
group, is: 



C^^.C^^lcpl'^'l-' (4.4.1.10) 



An outline of the proof, is as follows: 

Let two hexads 7^^^ and 7^^^ be given, such that their corresponding g's, are equal. 



Then there exists an automorphism matrix A (i.e. satisfying C^^^A^ = C^^-A^A^) 
connecting them, i.e. 7^^j = A^7^i''A^. But as it had been shown in the appendix of 



Iql , this kind of changes on 7^/3, can be seen to be induced by spatial diffeomorphisms. 
Thus, 3-dimensional Type II geometry, is uniquely characterized by some value of q. 

Although for full pure gravity, Kucharr |jl6|] has shown that there are not other 
first-class functions, homogeneous and linear in vr"'^, except Ha, imposing the extra 
symmetries (Type II), allows for such quantities to exist -as it will be shown. We are 
therefore, naturally led to seek the generators of these extra symmetries -which are 
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expected to chop off . Such quantities are, generally, called in the literature 

"Conditional Symmetries". 

The automorphism group for Type II, is described by the following 6 generators -in 
matrix notation and collective form: 



A(,)^ = ^ p \ (4.4.1.11) 

\ (T fl J 

with the property: 

C:.uK = + C:,\l (4.4.1.12) 

^From these matrices, we can construct the linear -in momenta- quantities: 

A(i) = Xli^plapT^"^ (4.4.1.13) 

Two of these, are the -ffa/s since C'^^-^^ correspond to the inner automorphism subgroup 
-designated by the x and y parameters, in A"^^^. The rest of them, are the generators of 
the outer automorphisms and are described by the matrices: 




=1 P \ (4.4.1.14) 

The corresponding -linear in momenta- quantities, are: 

= ^IdpIc^t:^^ (4.4.1.15) 
The algebra of these -seen as functions on the phase space, spanned by 7q/3 and vr'*'^-, is: 

[Ei, Ej} = CfjEx 

{Ei,H,} = -\XP^H0 (4.4.1.16) 
{Ei,H^} = -2(/€ + /i)7i? 



lYiom. the last of ( ^4.4.1.16| ), we conclude that the subgroup of i?/'s with the property 

K + /i = 0, i.e. the traceless generators, are first-class quantities; their time derivative 
vanishes. So let: 

Ei = {Ei: K + p = Q} (4.4.1.17) 
Then, the previous statement translates into the form: 

Ei = 0^ Ej = ci (4.4.1.18) 
the c/'s being arbitrary constants. 
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Now, these are -in principle- integrals of motion. Since, as we have earlier seen, 
Ej^s along with H^s, generate automorphisms, it is natural to promote the integrals of 
motion ( [4. 4.1.181) , to symmetries -by setting the c/'s zero. The action of the quantum 



version of these Ej^s on is taken to be 13 



Er^ = £?n.7„«^ = 



'ha = 



(/)/37ap^^^ U l^^^^j^g^^) (4.4.1.19) 



The Wheeler-DeWitt equation now, reads: 

5? ITT - st't^t + '^i^jnr + ^^ir - ^^ir - ^^^^ = o 4.4.1.20 

Note that: 

Vl = V' + ^^^R = (4.4.1.21) 

^"4(^-1) ^ ^ 

since we have a 2- dimensional, flat space, with contravariant metric: 

which is Lorentzian. This equation, can be easily solved by separation of variables; 
transforming to new coordinates u = and v = g7, we get the 2 independent equations: 

16u^A"{u) + 16uA'{u) - cA{u) =0 ^4 4 1 oq^ 

B"iv) + lB\v) - (^ + ^)B{v) = ^^•^•'•'^^ 



where c, is the separation constant. Equation ( |4.4.1.2(]| ), is of hyperbolic type and the 



resulting wave function will still not be square integrable. Besides that, the tracefuU 
generators of the outer automorphisms, are left inactive -due to the non vanishing CCR 
with Hq. 

These two facts, lead us to deduce that there must still exist a gauge symmetry, 
corresponding to some -would be, linear in momenta- first-class quantity. Our starting 
point in the pursuit of this, is the third of ( |4.4.1.16| ). It is clear that we need another 



quantity -also linear in momenta- with an analogous property; the trace of ti'^'^, is such 
an object. We thus define the following quantity: 

T = E,-(/€ + /i)7„^7r"^ (4.4.1.24) 

in the phase space -spanned by 7^/3 and tt^'^. It holds that: 

{T,Ho} = 

{T,Ha} = (4.4.1.25) 
{T,Ei} = 
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because of: 



{E,,7} = -2(«: + /i)7 

Again -as for Ej^s-, we see that since T, is first-class, we have that: 

T = ^ T = const = ct (4.4.1.27) 

another integral of motion. We therefore see, that T has all the necessary properties to 
be used in lieu of the tracefuU generator, as a symmetry requirement on \E'. In order to 
do that, we ought to set ct zero -exactly as we did with the c/'s, corresponding to i?/'s. 
The quantum version of T, is taken to be: 

Following, Dirac's theory, we require: 

dm d-m d^i d^i 

= A^7ap7^ + = + ^i){q^ - 1^) = (4.4.1.29) 

Equation ( ^.4.1.291 ), implies that ^1^(^,7) = "^{qj) and thus equation ( ^.4.1.201 ), finally, 
reduces to: 

4^2^" (w) + 4w^'(w) - 2w;^ = (4.4.1.30) 

where, for simplicity, w = qj. The solution to this equation, is: 

m = cMV2^) + C2Ko{^/2^) (4.4.1.31) 

where Jq is the modified Bessel function, of the first kind, and Kq is the modified Bessel 
function, of the second kind, both with zero argument. 

At first sight, it seems that although we have apparently exhausted the symmetries 
of the system, we have not yet been able to obtain a wave function on the space of the 
3-geometries, since \l/ depends on qj and not on q only. On the other hand, the fact that 
we have achieved a reduction to one degree of freedom, must somehow imply that the 
wave function found must be a function of the geometry. This puzzle finds its resolution 
as follows. Consider the quantity: 

Q = -27,.7r'"^ + ^^/^^^'-a7 7 lap7pa 40^^o,,7,^7 (4.4.1.32) 

q 

This can also be seen to be first-class, i.e. 

Q = O^Q = const = cn (4.4.1.33) 
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Moreover, it is a linear combination of T, -E/'s, and HaS, and thus = 0. Now it can 
be verified that Q, is nothing but: 

^ ^^ + -^) (4.4.1.34) 



So: 

-^gV3 = ^ = constant (4.4.1.35) 
Without any loss of generality, and since d is not an essential constant of the classical 



system (see [17] and reference [18] therein), we set "i? = 1. Therefore: 

m = cJo{V2q'/^) + C2Ko{V2q'/^) (4.4.1.36) 

where /q is the modified Bessel function, of the first kind, and Kq is the modified Bessel 
function, of the second kind, both with zero argument. 

As for the measure, it is commonly accepted that, there is not a unique solution. A 
natural choice, is to adopt the measure that makes the operator in ( [4.4.1.301) , hermitian 
-that is: 

/i(g)ocg-^ (4.4.1.37) 

It is easy to find combinations of ci and C2 so that the probability /i(g)|\E'p, be defined. 

Note that putting the constant associated with Q, equal to zero, amounts in restrict- 
ing to a subset of the classical solutions, since cq, is one of the two essential constants 
of Taub's solution. One could keep that constant, at the expense of arriving at a wave 
function with explicit time dependence, since then: 



7 



q-^/^Exp[ j cnN{t)dt] (4.4.1.38) 



We however, consider more appropriate to set that constant zero, thus arriving at a 
^ depending on q only, and decree its applicability to the entire space of the classical 
solutions. Anyway this is not such a blunder, since \E' is to give weight to all states, 
-being classical ones, or not. 
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The above treatise, remains almost unchanged for -(T. Chistodoulakis, G. O. Pa- 
padopoulos, to appear in Phys. Lett. B): 



4.2 The Class A Bianchi Type VI & VII Case 

The only difference between the various Bianchi Types, is that of the values of the 
structure constants. Of course, it must be emphasized that in the present work we study 
differential equations (i.e. the Einstein's Field Equations), which are defined locally -we 
are not concerned with global properties (i.e. topological aspects). 

It is well known that in 3 dimensions, the tensor admits a unique decomposition 
in terms of a contravariant symmetric tensor density of weight —1, m"^, and a covariant 
vector Va = ^C^^ as follows: 

= m'^'esf,^ + iyp6^ - iy^6^ (4.4.2.1) 
For the Class A (z/„ = 0) Bianchi Type VI {e = 1) and VII [e = —1), this matrix is: 

m"'^ = I -1 I (4.4.2.2) 




resulting in the following non vanishing structure constants 

C2^3 = e Cf, = 1 (4.4.2.3) 



In general, via the method of characteristics |jT2[, the quantum version of the three 
independent linear constraints can be used to reduce the dimension of the initial con- 
figuration space from 6 (jais) to 3 (combinations of 70/3), i-e. \E' = "if {q^ , q"^ , q^) [|^, 
where: 



la 



= C'^pAl^' (4-4.2.4) 



In analogy with the preceding case (where an outline of the proof is exhibited), the 
only G.C.T. (gauge) invariant quantities, which uniquely and irreducibly, characterize a 
3-dimensional geometry admitting the Class A Type VI and VII symmetry groups, are 
the quantities q^ and q^. 

In terms of the three q 's, one can define the following -according to |]13|^ induced 
"physical" metric, given by the relation: 



5g^g^ — IQq^q^ q^q^ q^q^ 



g'' = L^p,.^P^ = I qW qV 1 (4.4.2.5) 



q q'^ q q'^ -Sq'^q'^ 

Note that the first-class algebra satisfied by Hq, Ha, ensures that indeed, all components 
of g^^ are functions of the g*. 
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The automorphism group for the Class A Type VI and VII, is described by the 
following 4 generators -in matrix notation and collective form: 




Afn. = c a d \ (4.4.2.6) 



^From these matrices, we can construct the linear -in momenta- quantities: 

= XfDplapn''^ (4.4.2.7) 

Three of these are the H^s, since Cq^-j^ correspond to the inner automorphism subgroup 
-designated by the c, b and d parameters in A"^-j^. The remaining is the generator of the 
outer automorphisms and is given by the essentially unique matrix 

/ a \ 

e;^ = a (4.4.2.8) 
\ / 

The corresponding -linear in momenta- quantity is: 

E = e^TapTT''^ (4.4.2.9) 
It is straightforward to calculate the Poisson Brackets of E with Hq, Ha'. 

{E,Ha} = -lX^H(^ 

{E, Ho} = -2a^R = -2a^{q' + 2q^) [^A.Z.W) 

At this point, it is crucial to observe that we can construct a classical integral of 
motion, i.e. an extra gauge symmetry of the corresponding classical action: notice that 
the trace of the canonical momenta, 7^jy7r^'', has vanishing PB with Ha, E, and a -similar 
to E- PB with Hq equal to 2jR; thus, if we define: 

T^E- a-f^^n'"' (4.4.2.11) 

we can easily derive the following PB of T with Hq, Ha- 

{T,Ho} = 

{T,Hi} = -^Hi^O r4 4 2 12l 

{T,H,} = -IH2^0 [^A.z.LZ) 

{T,H,}^0 

The quantity T, is thus revealed to be first-class, and therefore an integral of motion 
(since the Hamiltonian, is a linear combination of the constraints): 

f ^ {T,H} !=^0^T ^ const = Ct (4.4.2.13) 
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The quantum version of T, is taken to be : 

(without any loss of generahty -see (|4.4.2.8D - we can safely suppress a, whenever conve- 
nient, by setting it to 1). 

Following the spirit of Dirac, we require: 

= e^7.p^ - ai^P^ = (q'ttj + q'jr-, - 7^) = Ct^ (4.4.2.15) 

The general solution [jl2l to the above equation has the form: 



^ = ^-t^T0(^^^g2) (4.4.2.16) 

being an arbitrary function in its arguments. 

Now, the number of our dynamical variables, is reduced from 3 (g*,z = 1,2,3) to 2, 
namely the combinations = q^ /q^ and w"^ = q^q^. So, we have a further reduction 
of the 3-dimensional configuration space spanned by the 3 q 's. Again, in terms of the 
w 's, the finally reduced "physical" -although singular- metric, is given by the following 
relation: 

i^dw^dw' /_i6 + 4(w;i)2 0\ 
^ =9''^^= n n (4.4.2.17) 



dq^ dqi V 

The singular character of this metric is not unexpected; its origin lies in the fact that 
I'^Pi&T^^T^^ = 0, where L"^^^ = 1(7"t^/35 ^^a<5^/37 _ 27"/^7T^) is the covariant superme- 
tric (inverse to Lajj-ys) and T^/j are the components of T, seen as vector field in the initial 
superspace spanned by 7^/3. Indeed, it is known that reducing to null surfaces entails all 
sort of peculiarities. 

So far, the degrees of freedom are two {w^, w"^). The vanishing of the s^^ component 
indicates that w"^ is not dynamical at the quantum level. This fact has its analogue at 
the classical level; indeed, consider the quantity: 

n=^^ = - L.,,.vr^- ^ (4.4.2.18) 

7g^ q^ 

In ( 4.4.2.1(^ ), the transition from velocity phase space to momentum phase space has 
been made using the Hamiltonian ( [4.4.1.3| ) -with the proper R. It is straightforward to 
verify that: 

n=-T- - (4.4.2.19) 
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Taking into account the weak vanishing of the hnear constraints, it is deduced that; 

2 2C 

n = -T = — - = 2Ct (4.4.2.20) 

a a 

if we set a = 1. 

Another way to show that VL is constant, is the following: 

{fi, H} = {^^^^^, H} = -^^{A^^^Hl + A^22Hl - 8e^i2H^H2) (4.4.2.21) 

Again, the weak vanishing of the linear constraints, ensures that Q is constant. 
Using (^4.4.2.181) , (|4.4.2.2q ) and the action, we have: 



_ Ciexp{2CT j N{t)dt} (4.4.2.22) 

Now, returning to the quantum domain, we observe that out of the three arguments 
of the wave function given in (|4.4.2.16| ), only q^/q^ is G.C.T. invariant -in the sense 



previously explained. This suggests that we must somehow eliminate 7, 7g^. To this 
end, we adopt the value zero for the classical constant Ct- This amounts to restricting 
to a 2-parameter subspace of the classical space of solutions, spanned by the 3 essential 



constants [0. This means that we base our quantum theory on this subspace and 
decree the wave function, to be applicable to all configurations (classical or not). The 
benefit of such an action, is twofold: 7"*^^ drops out, while at the same time w"^ = 7g^ 



is set equal to the constant Ci -see ( |4.4.2.22| ). These facts, along with the obligation 



that no derivatives with respect to w , are to enter the Wheeler-DeWitt equation -see 
( [4.4.2. 17[ )-, allow us to arrive at the following form for the wave function ( |4.4.2.16D : 



^ = e(4,Ci) (4.4.2.23) 

and of course: 

= q^q^ = const = Ci (4.4.2.24) 

Now, the final reduction of the configuration space is achieved. Our dynamic variable is 
the ratio q^/q"^, which is a combination of the only curvature invariants existing in Class 
A Bianchi Type VI and VII, and emerges as the only true quantum degree of freedom. 



Consequently, following the spirit of [|T3[ , we have to construct the quantum analogue 



of Hq as the conformal Laplacian, based on the non-singular part of the "physical" metric 



( [4.4.2.171) , i.e.: 



where: 



Ho = -lvl + w^{2 + w') (4.4.2.25) 



a.i{v^s^' 9^1} (4.4.2.26) 



/Sii 
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is the 1-dimensional Laplacian based on sn {s^^sn = 1). Note that in 1-dimension the 
conformal group is totally contained in the G.C.T. group, in the sense that any conformal 
transformation of the metric can not produce any change in the -trivial- geometry and 
is thus reachable by some G.C.T. Therefore, no extra term in needed in (|4.4.2.26 ), as it 



can also formally be seen by taking the limit d = 1, i? = in the general definition: 



= + -^^r^R = V (4.4.2.27) 



Thus, the Wheeler-DeWitt equation, reads: 



Cii2 + w) 



Ho^ = Viw^ - 4) d^Wiw^-A) 9^}^ - ^ = (4.4.2.28) 
where, for simplicity, we set = w and w'^ = Ci. 

Using the transformation w = 2cosh{z), the previous equation takes the form: 

- (Ci + cosh{z))^ = (4.4.2.29) 



The solutions to this family of equations are the Mathieu Modified Functions -see |]T 
and references therein, for an extended treatise, in various cases. 

As for the measure, it is commonly accepted that there is not a unique solution. A 
natural choice, is to adopt the measure that makes the operator in ( [4.4.2.28D hermitian, 
that is /i(w) = ^Ji_^ , or in the variable z, /i(z) = 1. However, the solutions to ( |4.4.2.29| ) 
can be seen to violently diverge for various values of 2 G [0, oo), which is the classically 
allowed region. If we wish to avoid this difficulty, we can abandon hermiticity, especially 
in view of the fact that we are interested in the zero eigenvalues of the operator, and thus 
does not make any harm to lose realness of the eigenvalues. If we adopt this attitude, we 
can find suitable measures, e.g. fi{z) = . The probability density p{z) = /i(2;)|\E'(z)p 
is now finite, enabling one to assign a number between and 1 to each Class A Type VI 
and VII geometry. 



Another feature of the reduction to ([4.4.2.29|) is that the final dynamical argument of 



\I' is the ratio q^/q^-, which is of degree zero in the scale factors, as seen from equation 



( [4.4.2.4|) . This consists a kind of build-in regularization with respect to the volume of 
the 3-space. Moreover, the solutions to this equation exhibit an increasingly oscillatory 
behavior, as Ci increases. This is most welcomed and expected in view of Ci being some 
kind of measure of the 3- volume, since it contains 7 -see (|4.4.2.24| ). 
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